DESIGN OF TRANSPIRATION COOLED THERMAL PROTECTION 

SYSTEMS 


E. Eugene Callens, Jr. 
Robert F. Yinet 


Mechanical Engineering 
Louisiana Tech University 
P.O. Box 10348 
Ruston LA 71272-0046 


March 1999 


Technical Report 


Contract No. 32-4136-58019 
and NASA(1997)-Stennis-17 
and NAS13-580 


Prepared for: 

National Aeronautics and Space Administration 
John C. Stennis Space Center 
Stennis Space Center, MS 39529-6000 


ABSTRACT 


This study explored three approaches for the utilization of transpiration cooling in thermal 
protection systems. One model uses an impermeable wall with boiling water heat transfer 
at the backface (Model I). A second model uses a permeable wall with a boiling water 
backface and additional heat transfer to the water vapor as it flows in channels toward the 
exposed surface (Model II). The third model also uses a permeable wall, but maintains a 
boiling condition at the exposed surface of the material (Model III). The governing 
equations for the models were developed in non-dimensional form and a comprehensive 
parametric investigation of the effects of the independent variables on the important 
dependent variables was performed. In addition, detailed analyses were performed for 
selected materials to evaluate the practical limitations of the results of the parametric study. 

It was found that for sufficiently thin walls using Model I, even matenals with low critical 
temperatures might be used in very high temperature environments. However, m many 
practical engineering applications, the thin walls do not provide adequate structure 
strength. Realization of the maximum theoretical gams for Model II requires sma 
diameter channels, larger material thickness, and higher environmental temperatures. Sma 
diameter channels typically have substantial problems with fouling which will limit their 
applicability. Above an environmental temperature of 3,600 K, it was found Modellll 
requires only moderately more water than the impermeable wall model (Model I). The 
maximum difference is approximately 19 percent at 5,000 K. Model m is the better design 
choice in this temperature range since a relative strong, ductile, and inexpensive materia 
such as steel can be used. In the range below 3,600 K, Model HI requires an increase in 
mass flow rate from 50 - 80 percent over niobium carbide. In this region, Model I using 
niobium carbide or some other suitable ceramic is a reasonable design alternative 

For the severe thermal environments such as the flame buckets of rocket engine test 
facilities, the boiling water exposed surface model (Model HI) is recommended in t e 
impingement regions. The impermeable wall model (Model I) with a layered combination 
of a ceramic, such as niobium carbide for the exposed material and steel or aluminum as 
the support material, is recommended in the peripheral areas. For thermal protection of 
personnel and equipment under conventional fire conditions, Model I with steel or 
aluminum as the wall material is recommended. The use of the impermeable wall with hea 
transfer to the vapor in the channels (Model II) should be limited to a restricted range o 
design requirements as its disadvantages outweigh its advantages for most materials an 
environmental conditions. 


l 



TABLE OF CONTENTS 


ABSTRACT 

LIST OF FIGURES 

LIST OF TABLES 

NOMENCLATURE 

Chapter 

1. INTRODUCTION 

Background 

Previous Work 

Approach of the Present Study 

2. MODEL I DEVELOPMENT AND APPLICATION ... 

Mathematical Model 

Parametric Study. 

3 . MODEL H DEVELOPMENT AND APPLICATION.. 

Mathematical Model 

Parametric Study 

4. MODEL m DEVELOPMENT AND APPLICATION 

Mathematical Model 

Comparison of Models I and HI 

5. CONCLUSIONS AND RECOMMENDATIONS 

Summary of General Approach 

Conclusions 

Recommendations 

REFERENCES 


Page 

i 

iv 

vii 

viii 


1 

1 

2 

3 

5 

5 

9 

23 

23 

28 


62 

63 

66 

66 

67 

69 

70 


li 



Page 


Appendix 

A PARAMETRIC STUDY OF THE AXIAL WALL 

AND VAPOR TEMPERATURE DISTRIBUTIONS 

B PARAMETRIC STUDY OF THE VARIATION OF THE RATIO 
OF THE VAPOR EXIT TEMPERATURE TO THE 
WALL SURFACE TEMPERATURE 

C PARAMETRIC STUDY OF THE VARIATION OF THE 

DIMENSIONLESS HEAT TRANSFER RATE, Q l2 

D. PARAMETRIC STUDY OF THE RATIO OF DIMENSIONLESS 

HEAT TRANSFER RATE FOR MODELS II AND I, Qu/ Qli 

E. PARAMETRIC STUDY OF THE RATIO OF DIMENSIONLESS 
MASS FLOW RATE PARAMETER, p, FOR MODELS II AND I. 

F. MAXIMUM ALLOWABLE ENVIRONMENTAL 

TEMPERATURE FOR MODEL I AS A FUNCTION OF 

MATERIAL THICKNESS AND CONVECTION HEAT 
TRANSFER COEFFICIENT (CROSS-PLOTS) 

G. PARAMETRIC STUDY OF THE RATIO OF DIMENSIONLESS 
MASS FLOW RATE PARAMETER, P, FOR MODELS 

n AND I - SELECTED MATERIALS 

H. PARAMETRIC STUDY OF MASS FLOW RATE PER UNIT 

AREA FOR MODELS I AND II 

I. COMPUTER PROGRAM DOCUMENTATION 


71 


80 


86 


92 


98 


104 


108 

136 

146 


iii 



LIST OF FIGURES 


1 . Schematic of Model I (Boiling Water Backface Model) 

2. Influence of Bi on T s /T b f as a function of T e /T b f 

3 . Influence of T e /T bf on T s /T bf as a function of Bi 

4. Influence of Bi on (T e /T bf ) ma x as a function of A 

5 . Influence of A on (T e /T bf )max as a function of Bi 

6. Schematic of Model II (Transpiration Model) 

7. Axial Wall and Vapor Temperature Distributions at the 

Baseline Conditions T e /T bf =6, [A=10 3 ], Bi=5, 7=2, a- 1.9742 

8. Variation of the Wall Surface Temperature and the 

Vapor Exit Temperature with Each Independent Variable 

9. Ratio of the Vapor Exit Temperature to the Wall Surface 

Temperature versus the Environmental Temperature 

10. Dimensionless Heat Transfer Rate, Q L2 , versus the 

Environmental Temperature 

1 1 . Ratio of Dimensionless Heat Transfer Rate for Models 

II and I, Q L2 / Qli, versus the Environmental Temperature 

12. Variation of Dimensionless Mass Flow Rate Parameter, P, 

for Models I and II with Each Independent Variable 

13. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models 

II and I versus the Environmental Temperature 

14. Maximum Allowable Environmental Temperature for Model I 
as a Function of Material Thickness and Convection 

Heat Transfer Coefficient 

1 5 . Comparison of Maximum Allowable Environmental Temperature 

for Model I as a Function of Material Thickness and 
Convection Heat Transfer Coefficient 


Page 

5 

11 

13 

16 

19 

23 

29 

30 
32 

. 34 

. 36 

. 39 

.. 41 


43 


47 


IV 



Page 


1 6. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Environmental Temperature and 

Material Thickness 

17. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models 13 
and I as a Function of Environmental Temperature and 

Channel Diameter 

1 8. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Environmental Temperature and 

Heat Transfer Coefficient 

19. Mass Flow Rate Per Unit Area for Models I and II as a Function 

of Environmental Temperature 

20. Schematic of Model HI (Surface Transpiration Cooling) 

2 1 . Comparison of Mass Flow Rate as a Function of Environmental 

Temperature for Model I and Model ID 

A1 . Parametric Study of the Axial Wall and Vapor Temperature Distributions 

B 1 . Parametric Study of the Variation of the Ratio of the V apor Exit 

Temperature to the Wall Surface Temperature 

C 1 . Parametric Study of the Variation of the Dimensionless 

Heat Transfer Rate, 

D1 . Parametric Study of the Ratio of Dimensionless Heat Transfer Rate 

for Models II and I, Qjj./ Qli 

El . Parametric Study of the Ratio of Dimensionless Mass Flow Rate 

Parameter, P, for Models II and I 

F 1 . Maximum Allowable Environmental Temperature for Model I 
as a Function of Material Thickness and Convection 
Heat Transfer Coefficient (Cross-Plots) 

Gl. Ratio of Dimensionless Mass Flow Rate Parameter, p, for 
Models II and I as a Function of Channel Diameter and 
Environmental Temperature 


49 


52 


55 

58 

62 

64 

72 

81 

87 

93 

99 


105 


109 


v 



Page 


G2. Ratio of Dimensionless Mass Flow Rate Parameter, p, for 
Models 13 and I as a Function of Channel Diameter and 
Material Thickness 

G3. Ratio of Dimensionless Mass Flow Rate Parameter, p, for 
Models II and I as a Function of Channel Diameter and 
Convection Heat Transfer Coefficient 

G4. Ratio of Dimensionless Mass Flow Rate Parameter, p, for 
Models 13 and I as a Function of Convection Heat Transfer 
Coefficient and Environmental Temperature 

G5. Ratio of Dimensionless Mass Flow Rate Parameter, p, for 
Models II and I as a Function of Convection Heat Transfer 
Coefficient and Material Thickness 

G6. Ratio of Dimensionless Mass Flow Rate Parameter, P, for 
Models II and I as a Function of Convection Heat Transfer 
Coefficient and Channel Diameter 

G7. Ratio of Dimensionless Mass Flow Rate Parameter, p, for 
Models II and I as a Function of Material Thickness and 
Environmental Temperature 

G8. Ratio of Dimensionless Mass Flow Rate Parameter, P, for 
Models 13 and I as a Function of Material Thickness and 
Channel Diameter 

G9. Ratio of Dimensionless Mass Flow Rate Parameter, p, for 
Models II and I as a Function of Material Thickness and 
Conduction Heat Transfer Coefficient 

HI . Mass Flow Rate Per Unit Area for Models I and II as a 
Function of Material Thickness 

H2. Mass Flow Rate Per Unit Area for Models I and II as a 
Function of Channel Diameter 

H3 . Mass Flow Rate Per Unit Area for Models I and II as a 
Function of Convection Heat Transfer Coefficient 


112 


115 


118 


121 


124 


127 


130 

133 

137 

140 

143 


vi 



LIST OF TABLES 


1 . Material Properties 


Page 

10 



NOMENCLATURE 


A Area 

Bi Biot Number, Eqn. (10) 

Ci Eqn. (49) 

C 2 Eqn. (48) 

C 3 Eqn. (47) 

C p Specific Heat at Constant Pressure 

h Heat Transfer Coefficient 

hf g Latent Heat of Vaporization 

H bf Eqn. (24) 

Hvi Eqn. (25) 

Hvii Eqn. (58) 

k Thermal Conductivity 

L Material Thickness 

m Mass Flow Rate 

P Perimeter of Flow Channel 

q Heat Transfer Rate 

Ql Eqn. (23) 

T Temperature 

x Direction Normal to Wall Surface 


Subscripts 

bf Backface Condition 
cond Conduction 
conv Convection 
e Environmental 
f Flame 

h Environmental Convection Condition 
rad Radiation 
s Surface 
trans Transpiration 
v Vapor 

w Wall 


Superscript 


Vlll 


Dimensionless Parameter 


Greek Symbols 


a Eqn. (32) 

Pi Eqn. (18) 

Pii Eqn. (59) 

y Eqn. (31) 

s Emissivity 

X 2 Eqn. (45) 

X 3 Eqn. (46) 

A Eqn. (9) 

a Stefan-Boltzman Constant 


CHAPTER 1 


INTRODUCTION 


Background 


Experience in the protection of reentry vehicle nosetips utilizing transpiration cooling can 
be applied to the development of transpiration cooled thermal protection systems for 
personnel and equipment. The high heat of vaporization of water permits protection against 
very large heat transfer rates with modest transpiration flow rates. The reentry vehicle 
designs that utilize this concept are able to protect metallic nosetips with an on-board water 
supply during the entire reentry trajectory for high reentry angles and, consequently, 
extremely high heat rates. If the system is designed to provide boiling at the surface, the 
temperature of the protected surface is limited to the boiling temperature of the water at the 
given pressure. It is easy to provide insulation for this relatively low temperature. The heat 
rates and exposure times that can be tolerated are limited only by the ability of the system 
to deliver the required flow rates. 


The most efficient model for using transpiration cooling is a thermal protection system 
design where the surface temperature is maintained at the maximum allowable operating 
temperature of the material. This permits the water to absorb the maximum heat for the 
thermal constraints as sensible heat in the liquid phase, latent heat in the phase change, and 
sensible heat in the vapor phase. Water has been found to be the most effective liquid for 
use as the coolant because of its relatively high heat capacity and heat of vaporization. In 
order to achieve this design goal, the outer protective layer must have a relatively high 
density of flow channels (high porosity) to provide sufficient area for the required heat 
transfer to the water within the material. 

Potential applications of transpiration cooled thermal protection systems at NASA-SSC 
include cooling of flame buckets; protection of fire fighters and their equipment; shields 
to permit rocket engine test and diagnostic equipment to be located near the high 
temperature exhaust plumes; and firewalls or shields to protect personnel and equipment 
from high temperature environments. 

The current and future planned use of existing rocket engine test facilities for ground 
testing of various size engines necessitates development of improved concepts for 
satisfying the cooling requirements for the flame deflectors. Research efforts are directed 
at providing a basic design methodology that includes materials and required flow rates for 
specific thermal environments. The basic design data can be utilized to predict system 
requirements for specific applications. One goal is to reduce the amount of water used, and 
the associated cost, during testing. 
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Previous Work 


Previous work at Louisiana Tech University in the use of transpiration cooling in thermal 
protection systems includes three unpublished senior design projects and one master s 
thesis (Ref. 1). The senior design projects include both experimental and theoretical studies 
of concepts based on transpiration cooling to provide protection for personnel and 
equipment in domestic fire environments. 

In the two experimental studies, the basic design consisted of flexible, transpiration cooled 
material that could be used as clothing for personnel or covering for equipment. Many 
material combinations were investigated with the most effective being an inner waterproof 
layer of aluminized Kevlar and an outer porous layer of Hexel fiberglass. Water was 
supplied at a controlled pressure to the interface between the materials. This design was 
tested in an oven at temperatures up to 550 2 C (1022 2 F). It was found that the maximum 
steady-state temperature for the water in the interface region at the maximum 
environmental temperature tested (550 2 C) was 46 2 C (1 15 2 F). It was concluded that this 
design would provide adequate protection for personnel and materials in domestic fire 
environments for extended periods of time. 

The other senior design project involved the design of a specific piece of fire-fighting 
equipment, namely, a transpiration cooled fire fighter’s glove. It was found in interviews 
with firemen that, using current gloves, their hands often become excessively hot and limit 
the amount of time that they can stay near certain fires. Major components in the design 
included two pressurized aluminum cylinders for water storage, tubing and insulation to 
transport the water to the glove and to distribute the water throughout the glove, and a 
transpiration medium to absorb water and allow the water to evaporate. It was found that 
in an environment with a 1300 e C (2372 Q F) flame temperature and a 700 -C (1292 F) 
environmental temperature, 1.8 kg (4 lbm) of water is required for a 30 minute exposure 
time. This is approximately 0.5 gallon or 2 liters of water which can easily be stored m a 
self-contained vessel and transported by a fireman. 

The objective of the master’s thesis by Gamblin (Ref. 1) was to investigate the potential use 
of thermal protection coatings, single tiles, and layered insulation systems to protect the 
walls of the flame buckets used in the testing of the Space Shuttle Main Engine (SSME) 
at NASA-SSC. The physical behavior of the hypothetical system was modeled by a plane 
wall boundary value problem with a convective frontface condition and a boiling water 
backface condition. The principal recommendations and conclusions follow: 

1 . Presently available thermal coatings were not recommended because of 
their lack of durability under severe thermal conditions. 

2. Metals were not recommended as the exposed tile material because of 
their relatively low melting temperatures or, in the case of tungsten, 
poor resistance to oxidation and brittleness. 
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3 . Refractory metals were not recommended as the exposed tile material 
because of inadequate service temperatures and/or poor resistance to 
oxidation. 

4. The recommended design was a two-layer composite with a ceramic 
material supported by a metal. The ceramics with the highest thermal 
performance were tantalum carbide, hafnium carbide, and niobium 
carbide, and the recommended metal support material was aluminum. 

5. The recommended design options were partial replacements of the 
existing transpiration cooled system with the designated two-layer 
composite system, resulting in reductions in cooling water flowrate of 
50 to 80 percent from the current value, depending on the magnitude of 
the replaced area in the flame buckets. 

A search of the Compendex, INSPEC, and CARL Uncover databases reveal very few 
studies of the application of transpiration cooling to the design of thermal protection 
systems. The most directly applicable of these investigations to the current study is an 
experimental study by Yano, Ochi, and Enya (Ref. 2). They designed, built, and tested a 
room with transpiration cooled exterior walls. The room was intended to serve as a safe 
haven for people during a fire. They showed that the interior walls of the room remained 
at the initial room temperature for a flame temperature of 800 e C (1472 2 F) for a test time 
of 25 minutes. All temperatures were at steady-state so it was expected that the 
transpiration walls would maintain the initial interior temperature indefinitely. 


Approach of the Present Study 


There are three approaches explored in the present study for the utilization of transpiration 
cooling for thermal protection. The first, entitled Model I, uses impermeable materials with 
water cooling of the inner surface of the system. Heat is absorbed by the water in the form 
of sensible heat to the liquid and latent heat to enable boiling, but the vapor is intentionally 
not superheated due to the reduction in the overall heat transfer coefficient associated with 
film-boiling. The temperature of the exposed surface is limited to the maximum operating 
temperature of the material. 

The second approach, entitled Model H, specifies a permeable outer layer which is exposed 
to the high heat transfer environment. In this model, a boiling condition is maintained at 
the backside of the protective material, and heat is also added to the water vapor as it flows 
in channels toward the exposed surface. The heat is absorbed by the water as sensible heat 
to the liquid, latent heat to enable phase change, and sensible heat to the vapor. As in 
Model I, the surface temperature of the exposed surface can be as high as the maximum 
operating temperature for the material. Model II results in a lower heat transfer rate to the 
surface when compared to a transpiration model where the surface is maintained at the 
boiling temperature due to the reduced temperature difference between the environment and 
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the surface. It also utilizes the maximum heat absorption capability of the water for the 
system temperature limits. 

The third approach, entitled Model HI, maintains a boiling condition at the surface of the 
exposed material. Heat is transferred to the water in the form of sensible heat to the liquid 
and latent heat to enable phase change. 

Model I is the same as Gamblin (Ref. 1) used in the flame bucket study. It is used here as 
a baseline model and applied to a range of environmental conditions from the domestic 
fires to the severe flame bucket thermal conditions. To this end, the model is developed 
in non-dimensional form and the influence of the dimensionless independent variables is 
studied parametrically. 

The comparative advantage of Model II is that it minimizes the water flow rate required to 
maintain the exposed surface at or below its maximum operating temperature. The 
disadvantages are the additional cost of fabrication of the permeable materials and the 
operational problems of maintaining unrestricted flow of the water through the channels 
under all operating scenarios. The advantages of Model I are lower fabrication costs and 
reduced operational problems. The disadvantages are higher required water flow rates for 
a design surface temperature and additional insulation requirements for the boiling area in 
some applications. Model III has the advantage of avoiding problems associated with high 
surface temperature such as reduced structural strength and oxidation. It has the 
disadvantage of having the largest water flow rates of the three approaches. This study 
provides quantitative results which can be used as a basis for selection of the most suitable 
approach for a given application. 

In both Model I and Model II the exposed surface temperature is maintained at or below 
the critical temperature for the material and the back face is maintained at a temperature 20 
K above the boiling temperature of water. This corresponds to a boiling heat transfer 
coefficient of approximately 30 kW/m 2 -K. 

Also, all three models are steady-state since preliminary calculations show that steady-state 
is reached in less than two minutes for most conditions and materials of interest in this 
study. The thermal protection systems being considered here are envisioned to be used for 
extended periods in excess of 10 minutes. 


4 



CHAPTER 2 


MODEL I DEVELOPMENT AND APPLICATION 
Mathematical Model 

In Model I, the material backface is maintained at the boiling condition. Heat is absorbed 
by the water in the form of sensible heat to the liquid and latent heat to enable boiling. 

The schematic of Model I is presented in Fig. 1. 



Figure 1. Schematic of Model I (Boiling Water Backface Model) 


The governing equation for the temperature distribution in the wall is the steady state heat 
conduction equation 


with boundary conditions as follows: 
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Applying the boundary conditions, the solution becomes 
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An alternate form is obtained by equating the heat transferred to the exposed surface to the 
heat transferred to the water 
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The independent variables in the above equations are A, Bi, x, and T e . The dependent 
variables are T(x), T s , Ql, and Pj. A parametric study is conducted in which the 
independent variables are varied over a range corresponding to typical materials and 
anticipated conditions. For the study, it is assumed that Th = T e which is a worst-case 
condition for a given maximum environmental temperature. It is also assumed that the 
backface temperature is maintained at 20 K above the boiling temperature. 


Parametric Study 

There were six materials selected for use as either exposed or support layers as shown in 
Table 1. This table presents the material thermal conductivity, melting point, and critical 
or maximum operating temperature which is defined for this study as 75 percent of the 
melting temperature. These materials were selected because they are a good representation 
of the available range in thermal conductivity and maximum operating temperatures. The 
ranges of the dimensionless variables contained in the parametric studies below correspond 
to this set of materials. Other materials are easily studied using the computer codes 
provided with this report. 
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Table 1. Material Properties 


Material 

Description 

k 

f BTU ^ 

k 

f W ) 

Tmp 

(°F) 

Tmp 

(K) 

Tcritical 

(K) 

[0.75 -T m p] 

[hr-ft-°F J 

(m-K J 

Aluminum* 

6061-T6 

90 

155.8 

1080 

855 

642 

Copper * 

pure 

227 

392.9 

1980 

135C 

1017 

Steel * 

AISIC1020 

27 

46.7 

2750 

1783 

1337 

Titanium f 


12.65 

21.9 

3056 

1953 

1465 

Tungsten f 


100.54 

174 

6128 

3660 

2745 

Ceramic ** 

Niobium Carbide 

17.34 

30 

6542 

3890 

2917 

* Ref. (3) 
t Ref. (4) 
** Ref. (i) 



This study covers the range of expected environmental (flame) temperatures from ambient 
up to 5 000 K. Heat transfer coefficients vary from the value of 0.03 kW/m -K for a 
standard industrial fire to 10 kW/m 2 -K for the maximum expected value in the flame 
buckets of the rocket test facilities. 

In Fig 2 the influence of Biot number on surface temperature is shown as a function of 
environmental temperature for four selected values of A, with the appropriate independent 
variable held constant. The non-dimensional parameter A may be thought of as a 
dimensionless material thickness and the Biot number as a dimensionless heat transfer 
coefficient. It is seen in Figs. 2a and 2b that at small material thickness and small heat 
transfer coefficients, the surface temperature is near the value of the backface temperature. 
However, as the heat transfer coefficient (Biot Number) increases, the surface temperature 
approaches the value of the environmental temperature. From Figs. 2c and 2d, it is 
observed that for larger material thickness, the surface temperature approaches the 
environmental temperature even for values of heat transfer coefficient near zero. 


The conclusion is that for sufficiently thin walls using a boiling water backface, even 
materials with low critical temperatures like aluminum can be used in very hig 
temperature environments. This concept is employed in a familiar Boy Scout 
demonstration where water is heated to boiling in a paper cup over hot coals. It is observed 
that the paper cup does not bum until all the water is boiled away because its surface 
temperature remains below the ignition temperature of the paper. It only works when the 
paper is sufficiently thin. 

For very high environmental temperatures in more practical engineenng applications, the 
thickness or more properly, thinness requirement, is often not achievable because of 
structural requirements. 
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Figure 2. Influence of Bi on T s /T bf as a function of T e /T bf 







These same effects are shown from a different perspective by cross-plotting the 
independent variables as shown in Fig. 3 where the influence of T e /T bf on T/T b f as a 
function of Bi for four different values of A is presented. As expected, it is seen in Fig. a 
that the surface temperature is equal to the backface temperature for all environmental 
temperatures when both the material thickness and heat transfer coefficient approach zero 
(A = 0 and Bi = 0). The curves tend to coalesce for large values of Bi corresponding to the 
asymptotic behavior observed in Fig. 2 for increasing Bi at constant T e /T bf . In Figs. 3c and 
d it is seen that the curves tend to coalesce for all values of Bi when the matenal thickness 
becomes sufficiently large. This corresponds to the observation that all of the constant 
T e /T bf curves tend toward horizontal lines in Figs. 2c and d. Furthermore, for combinations 
of large material thickness and large heat transfer coefficient, the surface temperature 
approaches the environmental temperature as shown by the nearly 45 degree slope on the 
high Bi lines in Fig. 3d. 


Figure 3. Influence of T e /T b f on T s /T b f as a function of Bi 
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(d) A =0.01 


i 



Bi 

20 

— 10 
5 


Te/T bf 

i 

| 


It is of interest to determine the maximum environmental temperature for which Model I 
can be used for specified values of material thickness and heat transfer coefficient. The 
maximum environmental temperature is obtained when the surface temperature is at its 
maximum or critical value. Fig. 4 shows the influence of heat transfer coefficient on 
(Te)max as a function of material thickness. The curves are presented in a semi-log format 
to facilitate differentiation of material thickness effects at low values of heat transfer 
coefficient. As expected from the parametric results above, the highest allowable values 
of environmental temperature are associated with the combination of low values of material 
thickness and heat transfer coefficient. The maximum allowable environmental 
temperatures approaches the critical surface temperature as the heat transfer coefficient 
increases for all values of material thickness. The value of the heat transfer coefficient for 
which the material thickness curves converge increases only modestly for higher critical 
surface temperatures. 
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Figure 4. Influence of Bi on (Te/Tbf)max as a function of A 
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In Fig. 5, the two independent variables are reversed to show the same effects from a 
different perspective. For the higher values of heat transfer coefficient, all of the curves 
coalesce toward a single curve as the material thickness increases. 

The conclusion from analyses of the results presented in Figs. 4 and 5 is that material 
selection for Model I is limited by the highest anticipated environmental temperature. 

Figure 5. Influence of A on (Te/T b f) max as a function of Bi 


(a) Aluminum 
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CHAPTER 3 


MODEL II DEVELOPMENT AND APPLICATION 


Mathematical Model 


In Model H, the material backface is maintained at the boiling condition as in Mode . 
However, channels in the wall conduct the vapor through the wall where it continues o 
absorb heat and approach the wall temperature in the limit 


The schematic of Model II is presented in Fig. 6. 



Figure 6. Schematic of Model II (Transpiration Model) 


The governing equations for the temperature distributions in the wall and the vapor flow 
are obtained from a heat balance over the wall and vapor flow elements. 
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~ h v P(T w -T v ) = 0 
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(28) 


x * T 

Normalize the governing equations with x* = — and T = — 
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Reduce the 2 nd order ODE system to a coupled 1 st order ODE system. 
Let 


* 

z 



( 35 ) 


Therefore 
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( 36 ) 
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dz 

dx 
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The boundary conditions follow: 

At x* = 0 : rj=l and T* = 1 


At x = 1 : 
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~dx 


= Ql 
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( 46 ) 
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An overall heat balance for the wall and coolant yields 
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Normalize the equation by multiplying through by 
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(56) 


(57) 


(58) 


Therefore 


_ Ql 

H hf + H v n 


(59) 


Model II requires the flow to be regulated so that all of the water is boiled at the backface 
and the vapor is superheated in the flow channels. This condition has the effect of coupling 
y and a through Pn so that 



l fg 


TbfCp v 


(60) 


The independent variables in the above equations are y, a, A, Bi, x, and T e . The dependent 
variables are T w (x), T v (x), T s , T ve , Ql, and Pn. 
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Parametric Study 


Similar to the Model I case, a parametric study is conducted in which the independent 
variables are varied over a range corresponding to typical materials and anticipated 
conditions. As before it is assumed that Th = T e . The baseline case for the parametric study 
is selected to be 

T e = 6 

A = 0.001 ( 61 ) 

Bi = 5 

y = 2 

Fig. 7 shows the axial temperature distribution through the wall and vapor for Model II and 
the linear axial temperature distribution through the wall for Model I, all at the baseline 
conditions. It is observed that the vapor flow reduces the temperature in the wall. The 
most efficient design that could be obtained for Model II is one in which the vapor exit 
temperature is equal to the wall surface temperature. This would correspond to a condition 
in which the full heat sink capacity of the vapor is realized for the existing temperature 
limits. The value of the vapor exit temperature relative to the surface temperature is a 
measure of efficiency rather than effectiveness since it is the actual exit temperature 
obtained compared to the maximum possible. In this sense, the difference between these 
temperatures at x = L (x* = 1) is indicative of the efficiency of Model D relative to Model 
I where a smaller difference corresponds to greater efficiency. The full parametric study 
of the wall and vapor temperature distributions is presented in the figures of Appendix A. 
The parametric study presents combinations of the following selected values of the 
independent variables: 

T e * = 2,4,6, 8 

A = 10' 5 , 10 4 , 10" 3 , 10' 2 ( 62 ) 

Bi = 0, 1,5,20 

y = 1,2, 3, 4 

It is observed that conditions for which the vapor exit temperature is near the limiting value 
of the wall surface temperature is associated with high values of y. Physically, y represents 
the effectiveness of the channels in transferring heat from the wall to the water vapor. 
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Figure 7. Axial Wall and Vapor Temperature Distributions at the Baseline 

Conditions 

T e /T b f=6, [A=10' 3 ], Bi=5, y=2, a=1.9742 



This effect is observed more directly in Fig. 8 where T s / T b f and T ve / Tbf are P| ot ^ ve ^ su ^ 
each of the independent variables with the remaining independent variables held constant 
at their baseline values. It is observed in Fig. 8d that the conditions associated with the 
higher efficiencies correspond to higher values of y as previously noted. However, when 
the plausible ranges for the independent variables are calculated, it is noted that the 
maximum value of a is approximately 5 which corresponds to a value of y of 
approximately 3.5 for these conditions. Therefore, the potential gam m efficiency 
suggested by Fig. 8d cannot be realized in practice for existing materials. 

Fig 8a shows that the efficiency of Model II decreases with increasing environmental 
temperature for these conditions. It is cautioned that observations for theseconihonsmay 
not be valid for other reference conditions and should not be generalized. Fig. 8 gg 
that the efficiency is nearly independent of A. In Fig. 8c, it is seen that lower values of Bio 
number are associated with higher efficiency which asymptotically approaches a constant 

value as Bi increases. 
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Figure 8. Variation of the Wall Surface Temperature and the Vapor Exit 
Temperature with Each Independent V ariable 


(a) Variation with Environmental Temperature 
A=T0' 3 , Bi=5, y=2 



(b) Variation with A 
Te/T bf =6, Bi=5, y=2 
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A more direct measure of the efficiency of Model II is found by plotting the ratio of the 
vapor exit temperature to the wall surface temperature versus the environmental 
temperature as shown in Fig. 9. In Fig. 9a the influence of A at the baseline conditions is 
illustrated. It is seen that efficiency decreases with increasing environmental temperature, 
and the effect of A is small. In Fig. 9b, higher efficiencies are observed to be associated 
with lower values of Biot number. The strong influence of higher values of y is seen in Fig. 
9c. These effects are further confirmed by the complete parametric variation of the 
independent variables at the baseline condition as presented in Appendix B. 


Figure 9. Ratio of the Vapor Exit Temperature to the Wall Surface Temperature 
versus the Environmental Temperature 


(a) Effect of A 
Bi=5, y=2 
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It is also of interest to study the heat transfer rate to the surface as a function of the various 
combinations of the independent variables. The dimensionless heat transfer rate for Model 
n, Ql 2 , is presented in Fig. 10 as a function of environmental temperature for the selected 
values of the remaining independent variables. Fig. 10a shows the expected increase in 
heat transfer rate with increasing environmental temperature and a moderate increase in 
heat transfer rate with increasing A. In Fig. 10b the heat transfer rate increases significantly 
with increasing Biot number at all values of environmental temperature. A very important 
result is shown in Fig. 10c in that the heat transfer rate becomes nearly independent of y for 
sufficiently high values of y at all values of environmental temperature. This probably is 
responsible for the previously observed increasing efficiencies with increasing y. The 
cross-plots for all the independent variables are shown in Appendix C. 


Figure 10. Dimensionless Heat Transfer Rate, Q L2 , versus the Environmental 

Temperature 


(a) Effect of A 
Bi=5, y=2 
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The heat transfer rate to the surface increases for Model II relative to Model I at the same 
environmental temperature due to the decreased wall surface temperature. This effect is 
shown in Fig. 11 as the ratio of the dimensionless heat transfer rate, Q 1 . 2 / Qli, versus 
environmental temperature for the same selected values of the remaining independent 
variables as presented in Fig. 10. In Fig. 1 la it is seen that the Model II heat transfer rate 
increases from approximately 5 to 40 percent over the Model I value as environmental 
temperature increases by a factor of 5. There is only a slight influence of A. In Fig. 10b 
the ratio of heat transfer rates increases moderately with increasing Biot number. The same 
near independence of y for sufficiently high values of y is observed in Fig 1 lc as was noted 
in Fig. 10c. The cross-plots of the ratio of heat transfer rates versus the other independent 
variables are shown in Appendix D. 

Figure 11. Ratio of Dimensionless Heat Transfer Rate for Models II and I, Q L 2 / 
Qli, versus the Environmental Temperature 


(a) Effect of A 
Bi=5, y=2 
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The most important of the output parameters for this study is the minimum mass flow rate 
required to maintain the backface temperature at its specified value. It is this flow rate that 
must be supplied by the thermal protection system for the given environmental conditions. 
In Fig. 12, the variation of the dimensionless parameter, P, with each independent variable 
is presented for Models I and E. The parameter, p, is proportional to the mass flow rate and 
is defined by Eqns. (1 8) and (56) as 


''-7 


B 

Pu A 



(18) 

K T b f 


L h fg 

(56) 

I k\v T b / 



It is observed in Fig. 12a that P increases with environmental temperature as expected, and 
the advantage of Model E over Model I (lower mass flow rate) also increases with 
increasing temperature. The parameter A, which is proportional to the material thickness, 
is seen to have only a small influence on p in Fig. 12b and the advantage of Model E 
remains approximately constant. It is observed in Fig. 12c that the value of the Bi has a 
significant influence on p for low values of Bi but above the value of approximately 5, p 
increases only slightly with increasing Bi. The channel heat transfer effectiveness 
parameter, y, is seen to have a major effect on the advantage of Model II over Model I in 
Fig. 12d. The ability of any particular design to achieve high values of this parameter is 
critical to the selection of Model E or Model I. 
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Figure 12. Variation of Dimensionless Mass Flow Rate Parameter, p, for Models I 

and II with Each Independent Variable 


(a) Variation with Environmental Temperature 
A=10' 3 , Bi=5, y=2 
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The potential advantage of Model II over Model I in terms of lower required mass flow rate 
is shown more directly in Fig. 13 where the ratio Ph/Pi is plotted versus the environmental 
temperature for selected values of the other independent variables. In Fig. 13a the potential 
advantage is seen to increase from approximately 10 percent at a temperature of twice the 
backface temperature to approximately 25 percent at a temperature of ten times the 
backface temperature. The parameter A has very little influence. A similar result is for the 
effect of Bi is observed in Fig. 13b. The crossover of the constant Bi curves is due to the 
fact that the convection heat transfer is more dominant at lower temperatures, and radiation 
is more dominant at higher temperatures. The strong influence of the parameter y is shown 
in Fig. 13c where a potential advantage of approximately 35 percent is observed at a 
temperature ratio of ten. Again this shows the importance of being able to achieve high 
values of y in a particular design. 

The cross-plots of the independent variables are presented in Appendix E. The conclusion 
is that the greater advantages of Model II over Model I is achieved with combinations of 
high y and high environmental temperatures with minor effects associated with A and Bi. 


Figure 13. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II and 
I versus the Environmental Temperature 


(a) Effect of A 
Bi=5, y=2 
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It is instructive to examine the actual values of maximum allowable environmental 
temperatures for the selected materials when the surface temperature is at the material s 
critical value. A baseline case was selected as follows: 

t. = 0.5 (r,_ +1) 

L = 12 mm ( 63 ) 

d = 3 mm 

h e = 1,000 W/m 2 -K 

In Fig. 14, two of these parameters are held at their baseline values while the other two are 
varied over their expected range. In Figs. 14a, b, c, and d it is seen that both material 
thickness and convection heat transfer coefficient have significant influence on the 
maximum allowable environmental temperature for aluminum, copper, steel, and titanium. 
The influence of the convection heat transfer coefficient on the maximum allowable 
temperature is observed to be small for both tungsten and niobium carbide in Figs. 14e and 
f. The influence of material thickness is only moderate for niobium carbide in Fig. 14f. The 
observed effects are further validated in the cross-plots of the independent variables 
presented in Appendix F. 


Figure 14. Maximum Allowable Environmental Temperature for Model I as a 
Function of Material Thickness and Convection Heat Transfer Coefficient 
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(b) Copper 
T cr it = 1016 K 



(c) Steel 
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(d) Titanium 
T crit = 1464 K 



(e) Tungsten 
T crit = 2745 K 
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(f) Niobium Carbide 
T cr it = 2917 K 

6000 
5000 
4000 
Te ’"“ 3000 

(K) 

2000 
1000 


0 0.005 0.01 0.015 0.02 0.025 

L (m) 



he 

(W/m K) 

— ♦ — 1 000 
— » — 5000 
10000 : 


A direct comparison of the maximum environmental temperatures is presented in Fig. 15a 
for a specified value of convection heat transfer coefficient of 1 ,000 W/m 2 -K. The ranking 
of these materials from highest to lowest is tungsten, copper, niobium carbide, steel, 
aluminum, and titanium. Copper ranks second in spite of its relatively low critical 
temperature because of its high thermal conductivity. 

Another direct comparison of the maximum environmental temperatures is presented in 
Fig. 15b for a specified value of material thickness of 0.012m (0.47 in). The ranking of 
these materials is the same as in Fig. 15a except that the curves for copper and niobium 
carbide cross at high values of the heat transfer coefficient as do the curves for aluminum 
and tungsten. 

When other important factors such as strength, weight, cost, and workability are considered, 
steel is an excellent choice among the metals for all but the most severe thermal conditions. 
The added weight, relatively high cost, and workability issues associated with tungsten 
limits its usability for most thermal environments. Since the ceramics are brittle and 
therefore not well suited for structural support, these results suggest that a layered 
combination of a ceramic such as niobium carbide for the exposed material and steel as the 
support material would be an effective design for severe thermal conditions. 


46 




Figure 15. Comparison of Maximum Allowable Environmental Temperature for 
Model I as a Function of Material Thickness and Convection Heat Transfer 

Coefficient 
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It is instructive to examine the actual values of the ratio (3n/pi for the selected materials of 
this study when the surface temperature is at the material’s critical value. The baseline case 
is the same as used in study of maximum allowable environmental temperature presented 
above, Eqn. (63). 

Figs. 16, 17, and 18 show this ratio versus environmental temperature for selected values 
of the three remaining variables for each of the materials. The advantage of Model II over 
Model I is less than 1 0 percent for aluminum, copper, and tungsten for all combinations of 
environmental temperature and the other variables. Steel approaches an advantage of 20 
percent and titanium 25 percent, with the greater advantages being at higher environmental 
temperature and material thickness. Niobium carbide also approaches an advantage of 
approximately 25 percent for higher temperature and larger thickness. 

A comprehensive parametric study involving all combinations of the four independent 
variables is presented in Appendix G. The conclusions are as follows: 

1 . Al umin um and copper show advantages of Model II over Model I of 
less than 10 percent for all combinations of the design variables. 

2. Tungsten shows advantages greater than 10 percent only for 
transpiration channels less than 2 mm in diameter. 

3. Steel shows a maximum advantage of 20 percent. 

4. Titanium and niobium carbide both show maximum advantages of 
approximately 30 percent. 

Realization of these maximum gains requires small diameter channels, larger material 
thickness, and higher environmental temperatures. The added cost and complexity of 
Model II relative to Model I must be considered in selecting the final design. 
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Figure 16. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II and 
I as a Function of Environmental Temperature and Material Thickness 


(a) Aluminum 
T cr it = 641 K 



(b) Copper 
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(c) Steel 
T cr i t = 1337 K 



(d) Titanium 
T cr it = 1464 K 
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Figure 17. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II and 
I as a Function of Environmental Temperature and Channel Diameter 
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Figure 18. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II and 
I as a Function of Environmental Temperature and Heat Transfer Coefficient 
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It is important to compare the dimensional mass flow rate per unit area as a function of 
environmental temperature for the selected materials as shown in Fig. 19. In these 
comparisons the surface temperature of each material is held at its critical value and the 
other independent variables are held at the baseline condition. 

It is observed that there is little difference between Models I and II for any of the materials. 
Also, for environmental temperatures up to 2,200 K, the lowest mass flow rates are for 
titanium, niobium carbide, and steel — in that order. Niobium carbide and steel give the 
lowest values for temperatures between 2,300 — 2,600 K. For temperatures between 2,600 
— 3,600 K, niobium carbide provides substantially lower flow rates than the other materials. 
Surprisingly, copper can handle temperatures between 3,600 — 4,200 K along with tungsten. 
From 4,200 - 5,000 K, tungsten is the only choice. However, as previously stated, tungsten 
has a severe oxidation problem that may preclude its selection. 

A more complete parametric study is presented in Appendix H. The conclusion is that 
reductions of mass flow rate up to 20 percent may be achieved for some materials at 
combinations of thick walls, small diameter channels, and high environmental 
temperatures. However, small diameter channels typically have substantial problems with 
fouling which will limit their applicability. 


Figure 19. Mass Flow Rate Per Unit Area for Models I and II as a Function of 

Environmental Temperature 
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(b) Copper 
T crit = 1016 K 
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CHAPTER 4 


MODEL HI DEVELOPMENT AND APPLICATION 


Mathematical Model 


In Model HI, there is sufficient flow in the channels to maintain the boiling condition at the 
exposed material surface. This avoids problems associated with a high surface temperature 
such as oxidation and loss of structural strength. However, there is a higher heat transfer 
rate from the environment due to the lower surface temperature and, therefore, more water 
flow is required. This chapter will explore the required increase in flow rate for Model EDI 
compared to Model I. The schematic of Model HI is presented in Fig. 20. 



The governing equation for the required mass flow rate to maintain boiling at the exposed 
surface is obtained from the heat balance between the radiation and convection from the 
environment to the water. 
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Comparison of Models I and III 

The only independent variable in Eqn. (68) is the environmental temperature. It is assumed 
that the surface temperature is maintained 20 K above the boiling which is 393 K for the 
standard sea level condition. A comparison of the mass flow rate versus environmental 
temperature for Model in compared to Model I for the selected materials of this study is 
shown in Fig. 2 1 . Model HI is the same for all materials as no material properties appear 
in Eqn. (68). For Model I, all materials are at the baseline condition, Eqn. (63), and the 
surface temperature of each material is its critical temperature. 

It is observed in Fig. 21a that above 3,600 K, Model KI requires only moderately more 
water than copper or tungsten. The maximum difference is approximately 19 percent at 
5,000 K. It is concluded that Model HI is the better design choice in this temperature range 
since a relative strong, ductile, and inexpensive material such as steel can be used. 

From Fig. 21b, it is observed that, in the range from 2,300 - 3,600 K, Model m requires 
an increase in mass flow rate from 50 — 80 percent over niobium carbide. In this region, 
Model I using niobium carbide or some other suitable ceramic may be a reasonable design 
alternative. 

Finally, Fig. 21c indicates that, in the range below 2,300 K, Model El requires nominally 
50 percent more mass flow than titanium or niobium carbide. Again, this is a sufficient 
difference to warrant consideration of Model I with either of these materials. 
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Figure 21. Comparison of Mass Flow Rate as a Function of Environmental 
Temperature for Model I and Model III 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 


Summary of General Approach 

There are three approaches explored in the present study for the utilization of transpiration 
cooling for thermal protection. They are summarized as follows: 

• Model I has an impermeable wall with a boiling water backface. The 
temperature of the exposed surface is limited to the maximum operating 
temperature of the material. 

• Model II has a permeable wall with a boiling water backface. However, 
in this model, heat is also added to the water vapor as it flows in 
channels toward the exposed surface. As in Model I, the temperature 
of the exposed surface can be as high as the maximum operating 
temperature of the material. 

• Model III has a permeable wall with a boiling condition at the surface 
of the exposed material. The exposed surface is therefore maintained 
at a temperature 20 K above the boiling temperature of the water. 

The comparative advantage of Model II is that it minimizes the water flow rate required to 
maintain the exposed surface at or below its maximum operating temperature. The 
disadvantages are the additional cost of fabrication of the permeable materials, the 
operational problems of maintaining unrestricted flow of the water through the channels 
under all operating scenarios, and problems associated with high surface temperature. 

The advantages of Model I are lower fabrication costs and reduced operational problems. 
The disadvantages are higher required water flow rates for a design surface temperature, 
additional insulation requirements for the boiling area in some applications, and problems 
associated with high surface temperature. 

Model HI has the advantage of avoiding problems associated with high surface temperature 
such as reduced structural strength and oxidation. It has the disadvantage of having the 
largest water flow rates of the three approaches. 

This study provides quantitative results, which can be used as a basis for selection of the 
most suitable approach for a given application. A computer program is provided with this 
report that allows the user to calculate the surface temperature and the water mass rate for 
any selected material properties and thermal conditions. The description of the program 
and user instructions are presented in Appendix I. 
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Conclusions 


In this study, the governing equations for the models were developed in non-dimensional 
form and a comprehensive parametric investigation of the effects of the independent 
variables on the important dependent variables was performed. In addition, detailed 
analyses were performed for selected materials to evaluate the practical limitations of the 
results of the parametric study. The conclusions are classified by model type as follows: 

Model I: 

• For sufficiently thin walls using a boiling water backface, even materials 
with low critical temperatures like aluminum can be used in very high 
temperature environments. However, for very high environmental 
temperatures in more practical engineering applications, the thickness 
or more properly, thinness requirement, is often not achievable because 
of structural requirements. 

• Material selection for Model I is limited by the highest anticipated 
environmental temperature. The ranking of these materials from highest 
to lowest is tungsten up to 5,000 K, copper up to 4,000 K, niobium 
carbide up to 3,800 K, steel up to 2,700 K, aluminum up to 2,300 K, 
and titanium up to 2,200 K. Copper ranks second in spite of its 
relatively low critical temperature because of its high thermal 
conductivity. In practice, the use of these materials at these 
temperatures may be limited by unacceptable oxidation. 

Model II: 

From the non-dimensional parametric analyses, the advantages of Model El over Model I 
in terms of reduction of required mass flow rate are summarized as follows: 

• Aluminum and copper show advantages of Model II over Model I of 
less than 10 percent for all combinations of the design variables. 

• Tungsten shows advantages greater than 10 percent only for 
transpiration channels less than 2 mm in diameter. 

• Steel shows a maximum advantage of 20 percent. 

• Titanium and niobium carbide both show maximum advantages of 
approximately 30 percent. 

When the analyses were conducted for the selected materials over the entire range of 
anticipated environmental conditions, the following conclusions were reached: 
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• Realization of the above maximum gains requires small diameter 
channels, larger material thickness, and higher environmental 
temperatures. However, small diameter channels typically have 
substantial problems with fouling which will limit their applicability. 

The added cost and complexity of Model II relative to Model I must be 
considered in selecting the final design. 

• When other important factors such as strength, weight, cost, and 
workability are considered, steel is an excellent choice among the 
metals for all but the most severe thermal conditions. 

• The added weight, relatively high cost, and workability issues associated 
with tungsten limits its usability for most thermal environments. 

• Since the ceramics are brittle and therefore not well suited for structural 
support, these results suggest that a layered combination of a ceramic 
such as niobium carbide for the exposed material and steel or aluminum 
as the support material would be an effective design for severe thermal 
conditions. 

Model III: 

In Model ID, there is sufficient flow in the channels to maintain the boiling condition at the 
exposed material surface. This avoids problems associated with a high surface temperature 
such as oxidation and loss of structural strength. However, there is a higher heat transfer 
rate from the environment due to the lower surface temperature and, therefore, more water 
flow is required. This study established these water flow requirements as follows: 

• Above 3,600 K, Model EDI requires only moderately more water than 
Model I. The maximum difference is approximately 19 percent at 5,000 
K. It is concluded that Model m is the better design choice in this 
temperature range since a relative strong, ductile, and inexpensive 
material such as steel can be used. 

• In the range from 2,300 — 3,600 K, Model HI requires an increase in 
mass flow rate from 50 — 80 percent over niobium carbide. In this 
region, Model I using niobium carbide or some other suitable ceramic 
may be a reasonable design alternative. 

• Below 2,300 K, Model III requires nominally 50 percent more mass 
flow than titanium or niobium carbide. Again, this is a sufficient 
difference to warrant consideration of Model I with either of these 
materials. 
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Models II and in have an additional potential advantage that was not evaluated in this 
study. When the water vapor is injected into the high temperature thermal environment, 
it will have a cooling effect. Specifically, this is expected to reduce the convection heat 
transfer coefficient. For large injection rates, this effect could be si gni ficant. 


Recommendations 


Based on the conclusions of this study, the following recommendations are made for 
application of these models to the NASA-SSC thermal protection problems: 

• For the severe thermal environment of the flame buckets. Model in 
should be used in the impingement regions. Model I with a layered 
combination of a ceramic, such as niobium carbide for the exposed 
material and steel or aluminum as the support material, should be used 
in the peripheral areas. 

• For thermal protection of personnel and equipment under conventional 
fire conditions, Model I should be used with steel or aluminum as the 
wall material. 

The use of Model II should be limited to a restricted range of design requirements as its 
disadvantages outweigh its advantages for most materials and environmental conditions. 
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APPENDIX A 


PARAMETRIC STUDY OF THE AXIAL WALL AND VAPOR TEMPERATURE 

DISTRIBUTIONS 
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(g) Te/T bf =6, [A=ir 2 ], Bi=5, y=2, a=L7816 
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APPENDIX B 


PARAMETRIC STUDY OF THE VARIATION OF THE RATIO OF THE VAPOR 
EXIT TEMPERATURE TO THE WALL SURFACE TEMPERATURE 


80 




Figure Bl. Parametric Study of the Variation of the Ratio of the Vapor Exit 
Temperature to the Wall Surface Temperature 


(a) Influence of Environmental Temperature and A 
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APPENDIX C 


PARAMETRIC STUDY OF THE VARIATION OF THE DIMENSIONLESS HEAT 

TRANSFER RATE, Q L2 
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Figure Cl. Parametric Study of the Variation of the Dimensionless Heat Transfer 

Rate, Ql 2 


(a) Influence of Environmental Temperature and A 
Bi=5, y=2 
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APPENDIX D 


PARAMETRIC STUDY OF THE RATIO OF DIMENSIONLESS HEAT TRANSFER 
RATE FOR MODELS H AND I, Q L2 / Qli 
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Figure Dl. Parametric Study of the Ratio of Dimensionless Heat Transfer Rate for 

Models II and I, Q L2 / Qlx 
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(g) Influence of Environmental Temperature and y 
A=10' 3 , Bi=5 
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APPENDIX E 


PARAMETRIC STUDY OF THE RATIO OF DIMENSIONLESS MASS FLOW 
RATE PARAMETER, p, FOR MODELS H AND I 
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Figure El. Parametric Study of the Ratio of Dimensionless Mass Flow Rate 
Parameter, P, for Models II and I 

(a) Influence of Environmental Temperature and A 
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APPENDIX F 


MAXIMUM ALLOWABLE ENVIRONMENTAL TEMPERATURE FOR MODEL I 
AS A FUNCTION OF MATERIAL THICKNESS AND CONVECTION HEAT 
TRANSFER COEFFICIENT (CROSS-PLOTS) 
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(e) Tungsten 
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APPENDIX G 


PARAMETRIC STUDY OF THE RATIO OF DIMENSIONLESS MASS FLOW 
RATE PARAMETER, (3, FOR MODELS D AND I - SELECTED MATERIALS 
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Figure Gl. Ratio of Dimensionless Mass Flow Rate Parameter, J3, for Models II 
and I as a Function of Channel Diameter and Environmental Temperature 
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Figure G2. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Channel Diameter and Material Thickness 
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Figure G3. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Channel Diameter and Convection Heat Transfer 

Coefficient 
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Figure G4. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Convection Heat Transfer Coefficient and Environmental 

Temperature 
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Figure G6. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Convection Heat Transfer Coefficient and Channel 

Diameter 
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Figure G7. Ratio of Dimensionless Mass Flow Rate Parameter, J3, for Models II 
and I as a Function of Material Thickness and Environmental Temperature 
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Figure G8. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Material Thickness and Channel Diameter 
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Figure G9. Ratio of Dimensionless Mass Flow Rate Parameter, p, for Models II 
and I as a Function of Material Thickness and Conduction Heat Transfer 

Coefficient 
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PARAMETRIC STUDY OF MASS FLOW RATE PER UNIT AREA FOR MODELS I 

AND IT 
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Figure HI. Mass Flow Rate Per Unit Area for Models I and II as a Function of 

Material Thickness 
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Figure H3. Mass Flow Rate Per Unit Area for Models I and II as a Function of 
Convection Heat Transfer Coefficient 
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APPENDIX I 


COMPUTER PROGRAM DOCUMENTATION 
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The analytical models have been solved using a computer routine. The source code was 
developed and programmed in FORTRAN, and is divided into three files: MAIN. FOR, 
DEC_VAR . INC, and INPUT . DAT. The MAIN . FOR file is the program driver. The 
DEC_VAR . INC file contains all of the variable declarations. The INPUT . DAT file is the 
descriptive input deck for the models, and can be adjusted using any text editor. 

The MAIN. FOR and DEC_VAR.INC files have been compiled into an executable file, 
DTCTPS . EXE. This file can be executed either from a Windows environment or from MS- 
DOS. The DTCTPS . EXE file must be located in the same directory as the INPUT . DAT file. 
The program output will be to the screen, and contains both the relevant characteristic input 
data and the respective Model results. A data output file can be created by using the following 
command from an MS-DOS prompt: DTCTPS > OUTPUT. OUT. 

In order to solve Model I and II, a secant iteration scheme was used. Each secant iteration 
requires two initial guesses that can either be bounded or unbounded, and a numerical 
convergence criterion. All of these parameters are accessible in the INPUT . DAT file. Model 
II requires extra care in setting the initial guesses, as it has a highly non-linear iterative 
solution procedure. The default initial values of xL_2, xR_2, xL_3, and xR_3 should 
handle most all input conditions, however, they may need to be adjusted for the Model II 
solution to converge. The need for such adjustments may be indicated by infinite run-times 
and dubious results. In these situations, either decrease the convergence criteria or define 
more accurate initial guesses. 

All files are presented as follows: 
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MAIN. FOR 


PROGRAM DTCTPS 


* =z = = = = = = = = = = = = = = = = = = = = = = = z= = — ====================================== =====* 

*====—=================================================================* 

* Design of Transpiration Cooled Thermal Protection Systems * 

* NASA Stennis Space Center * 

* - * 

* E. Eugene Callens, Jr. & Robert F. Vinet * 

* Louisiana Tech University * 

* Ruston, Louisiana * 

* - * 

* March 1999 * 


* March 1999 



c Variables 

IMPLICIT NONE 
INCLUDE ' DEC„VAR . INC 1 


c Equations 

DOUBLE PRECISION EQ1, EQ2 
EXTERNAL EQ1, EQ2 


c 


Input 


NAMELIST 

NAMELIST 

NAMELIST 

NAMELIST 

OPEN 

(8, 

READ 

(8, 

READ 

(8, 

READ 

(8, 

READ 

(8, 


/MISC/ sigma, epsilon, hfg, Cp_liquid, Cp_vapor, 
Tinitial, Tbf, Te_r, Te_jh., hbar_e 
/MATERIAL/ MatName, kw, L 
/M0DEL2_0NLY/ hbar_v, P, Aw 

/ ITER_PAR/ EPS 1 , xL_l, xR_l , EPS_2, xL_2 , xR„2 , 

EPS„3, xL_3, xR_3 

FILE = ' input . dat 1 , FORM = 'FORMATTED 1 ) 

MISC) 

MATERIAL) 

M0DEL2_0NLY ) 

ITER_PAR) 


Tboil, 


************************************************************************ 

* STEADY STATE MODEL I * 

************************************************************************ 

c Call to the secant iteration to find value of Ts_l 
Falpha_l = 0 . 0D0 

CALL SITER (EQ1 , xL_l , xRJ. , Falpha_l, EPS_1, RESULT„1) 

Ts 1 = RESULT_1 

mDOTperA„l = ( kw* ( (Ts_l -Tbf ) /L) ) * 

$ { 1/ (hfg + Cp„liquid* (Tboil -Tinitial) 

$ + Cp„vapor* (Tbf -Tboil) ) ) 

************************************************************************ 

* STEADY STATE MODEL II * 

************************************************************************ 


c Call to the secant iteration to find value of alpha (& Beta_2) 
Falpha_2 = 0 . 0D0 

CALL SITER (EQ2 , xL_2 , xR_2 , Falpha_2, EPS_2 , RESULT_2) 
alpha = RESULT_2 
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mD0TperA_2 = (Beta_2 * (kw * Tbf ) ) / (L * hfg) 

************************************************************************ 

* STEADY STATE MODEL III * 

************************************************************************ 

c Set Ts at 20 K above Boiling Temp, at Stnd. Atm. Conditions 
Ts_3 = 393 

mDOTperA_3 = ( ( (sigma*epsilon) * ( (Te_r**4 ) - (Ts_3 **4 ) ) ) + ( hbar„e * 

$ (Te_h-Ts__3 ) ) ) * (1/ (hf g+Cp_liquid* (Tboil -Tinitial) ) ) 


* OUTPUT * 


WRITE (*, *) 

1 * 

1 

/ 

Mat Name, **■ 





WRITE (*, *) 

1 * 







WRITE (*,25) 

I * 

t 

1 kw ~ ‘ / 

kw. 

' W/m*K' , 

i * i 


WRITE (*,27) 

t * 

/ 

•i* = ■ , 

L, 

'm' , 

! * 1 


WRITE (*,26) 

1 * 

/ 

• Te_r = 1 , 

Te_r 

'K' , 

1 * 1 


WRITE (*,26) 

> * 

/ 

1 Te__h = 1 , 

Te_h 

'K' , 

1 * 1 


WRITE (*,26) 

1 * 

/ 

’Tbf = 1 , 

Tbf, 

'K' , 

1 * 1 


WRITE (*,28) 

» * 

/ 

1 hbar_e = 1 , 

hbar_ 

_e, 'W/m A 2*K' , 

1 * 1 


WRITE (*,28) 

1 * 

/ 

'hbar v = 1 , 

hbar_ 

_v, 'W/m A 2*K' , 

1 * 1 


WRITE (*,27) 

1 * 

/ 

, P = *, 

p. 

'm* , 

1 * 1 


WRITE (*,29) 

1 * 

t 

'Aw = 1 , 

Aw, 

1 m A 2 ' , 

! * 1 


WRITE (*,*) 1 

1 * _ . 







1 * \ 

t 

' MODEL I ' , 

i * i 




WRITE (*,32) 

WRITE ( * , 33 ) 

1 * 1 

/ 

*Ts_I 

- ' 

! , Ts_l , 

' K ' , 

« * i 

WRITE (*,34) 

1 * 1 

f 

, mdot/A„I 

= i 

, mDOTperA_l , 

' kg/m A 2*s ' , 

i * i 

WRITE (*, 32) 

1 * 1 

t 

'MODEL II ' , 

i * i 




WRITE (*,33) 

1 * 1 

/ 

' Ts_II 

— i 

, Ts_2 , 

1 K 1 , 

i * i 

WRITE (*,33) 

1 * 1 

/ 

' Tve_II 

= i 

, Tve_2 , ' 

'K', 

» * i 

WRITE (*,34) 

1 * I 

/ 

'mdot/A_II 

— i 

, mDOTperA„2, 

1 kg/m A 2*s 1 , 

i * » 

WRITE (*,32) 

1 * I 

r 

'MODEL III', 

i * i 




WRITE (* , 33)' 

1 * > 

i 

' Ts_III 

= • 

, Ts_3 , 

1 K 1 , 

i * i 

WRITE (*,34) 

1 * t 

/ 

1 mdot/A_III = ' 

, mDOTperA„3 , 

'kg/mA2*s ' , 

i * i 




25 

FORMAT 

(IX, 

A, IX, 

A, 

F6.1, IX, A, 29X, A) 

26 

FORMAT 

(IX, 

A, IX, 

A, 

F6.1, IX, A, 33X, A) 

21 

FORMAT 

(IX, 

A, IX, 

A, 

F6.4, IX, A, 33X, A) 

28 

FORMAT 

(IX, 

A, IX, 

A, 

F7.1, IX, A, 26X, A) 

29 

FORMAT 

(IX, 

A, IX, 

A, 

F10.8, IX, A, 27X, A) 

32 

FORMAT 

(IX, 

A, IX, 

A, 

4 IX, A) 

33 

FORMAT 

(IX, 

A, 14X, 

A, 

F6.1, IX, A, 16X, A) 

34 

FORMAT 

(IX, 

A, 14X, 

A, 

F8.4, IX, A, 7X, A) 


STOP 

END 


************************************************************************ 

* EQUATION 1: Steady State Model I of Ts equation (set equal to zero)* 

************************************************************************ 

DOUBLE PRECISION FUNCTION EQ1 (Ts_l) 

IMPLICIT NONE 
INCLUDE 1 DEC_VAR . INC 1 
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EQ1 = ( ( (sigma*epsilon*L) /kw) * ( (Te„r**4) - (Ts_l**4) ) ) + 
$ ( (hbar_e*L) /kw) * (Te_h - Ts_l) + Tbf - Ts__l 


END 

************************************************************************ 

* EQUATION 2: Steady State Model II (iteration on alpha) * 

************************************************************************ 

DOUBLE PRECISION FUNCTION EQ2 (alpha_iter) 

IMPLICIT NONE 
INCLUDE 1 DEC_VAR . INC 1 

c Equations 

DOUBLE PRECISION EQ3 
EXTERNAL EQ3 

alpha = alpha_iter 

C Calculate values 

gamma = (hbar_v * P * (L**2) ) / (kw * Aw) 

lambda2 = (-alpha + SQRT ( (alpha**2) +4*gamma ) ) / 2 

lambda3 = (-alpha - SQRT ( (alpha**2) +4*gamma ) ) / 2 

c Call to the secant iteration to find value of QL2 
Falpha„3 = 0.0D0 

CALL SITER (EQ3, xL„3 , xR„3 # Falpha„3 , EPS_3 , RESULTJ) 

QL2 = RESULT_3 

Tve_2 = (Cl + { (alpha/ (alpha+lambda2 ) ) * C2*EXP (lambda2 ) ) 

$ + ( (alpha/ (alpha+lambda3) ) * C3*EXP (lambda3) )) * Tbf 

Hbf = (hfg + Cp_l i quid* Tbf * ( (Tboil/Tbf ) - (Tini tial/Tbf ) ) ) /hf g 
Hv_2 = (Cp__vapor*Tbf * ( (Tve_2/Tbf ) - (Tboil/Tbf) ) ) /hfg 

Beta„2 = QL2/ (Hbf+Hv_2) 

EQ2 = ((gamma/alpha) * (hf g/ (Cp_vapor*Tbf ) ) ) - Beta_2 

END 

************************************************************************ 

* EQUATION 3: Steady State Model II (iteration on QL2) * 

************************************************************************ 

DOUBLE PRECISION FUNCTION EQ3 (QL2_iter) 

IMPLICIT NONE 
INCLUDE ' DEC_VAR . INC 1 

QL2 = QL2_iter 

c Calculate values 

C3 = ( (alpha+lambda3) *QL2) 

$ / (lambda3* (EXP ( lambda3 ) * (alpha+lambda3 ) 

$ - EXP (lambda2) * (alpha+lambda2 ) ) ) 

C2 = -C3* ( (lambda3* (alpha+lambda2) ) 

$ / (lambda2* (alpha+lambda3) ) ) 

Cl = 1 - C2 - C3 

Ts_2 = (Cl + C2*EXP(lambda2) + C3*EXP (Iainbda3) ) * Tbf 

EQ3 = {(sigma * epsilon * L * (Tbf **3) ) /kw) * ( (Te_r/Tbf ) **4 
$ - (Ts_2/Tbf) **4 ) + ( (hbar_e * L) /kw) * ( (Te_h/Tbf ) 

$ - (Ts_2/Tbf ) ) - QL2 
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END 


* SUBROUTINES * 

* * 

* Robert F. Vinet 23 March 1998 * 

* * 

* 

* 


* 

The Secant Iteration Function * 

* * 

************************************************************************ 


* This scheme uses the Secant Method Iteration for non-linear 

* equations as presented in Hoffman's book pg. 101. [Numerical 

* Methods for Engineers and Scientists, McGraw-Hill, 1992]; 

* The two initial guesses can either be bounded or unbounded. 

* 


* . 

* 

* 

* 

* 

* XL & XR 

* F ALPHA 

* 

* 

* 

* 

* 

* 


* 

* 

* 

* 

* 

- * 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 


EQ1 

EPS 


XC 


FL, FR, FC 


EQn : adjust what equation is sent into the iterator 

in the main program 

: defines the convergence criterion for the 

independent variable 
: two initial guesses required 

: sets value of f (x) to be solved; typically let 

FALPHA = 0 


intermediate value of the independent variable; 
a function of the type of iterator used 
name of the function to be iterated 
function evaluated at its respected ind. variable 


* * 
************************************************************************ 
SUBROUTINE SITER (F, XL, XR, FALPHA, EPS, RESULT) 

IMPLICIT NONE 
INTEGER COUNT 

DOUBLE PRECISION F, FALPHA, XL, XR, EPS, XC, FC, FL, FR, RESULT 
COUNT = 0 
FR = F (XR) 

FL = F (XL) 

XC = ( (FALPHA* (XL- XR) ) + (FL*XR) - (FR*XL) ) / (FL - FR) 

FC = F (XC) 

10 IF (ABS(XL-XR) . LE . (2.0*EPS)) GO TO 100 

. GT . 0 ) THEN 


IF ( (FC* 

FL) 

XR 

- 

XL 

FR 

= 

FL 

XL 

= 

XC 

FL 

= 

FC 

ELSE 



XR 

= 

XL 

FR 

= 

FL 

XL 

= 

XC 

FL 

= 

FC 


ENDIF 

XC = ( (FALPHA* (XL- XR) ) + 
FC = F (XC) 

COUNT = COUNT + 1 
GOTO 10 


(FL*XR) - (FR*XL) ) / (FL-FR) 
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C100 WRITE (7,*) 'Number of Iterations to Converge = COUNT 
100 RESULT = ( (FALPHA* (XL-XR) ) + (FL*XR) - (FR*XL) ) / (FL-FR) 
RETURN 
END 
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DEC_VAR . INC 


c Input Variables 


DOUBLE PRECISION sigma, epsilon, hfg, Cp liquid, Cp_vapor, Tboil 

DOUBLE PRECISION Tinitial, Tbf, Te_r, Te_h, hbar_e, kw, L, hbar_v 
DOUBLE PRECISION P, Aw 
CHARACTER* 50 Mat Name 

c Output Variables 


DOUBLE PRECISION Ts_l , Ts_2 , Ts_3 , mDOTperA_l , mDOTperA_2 
DOUBLE PRECISION mDOTperA_3 , Tve_2 

c Intermediate & Working Variables 


DOUBLE PRECISION EPS_1 , xL_l , xR_l , Falpha_l, RESULT_1 

DOUBLE PRECISION EPS_2 , xL_2 , xR_2 , Falpha_2, RESULT_2 

DOUBLE PRECISION EPS_3 , xL_3 , xR_3 , Falpha_3 , RESULT_3 

DOUBLE PRECISION gamma, lambda2, lambda3, alpha, QL2 , Cl, C2, C3 

DOUBLE PRECISION Hbf , Hv_2 , Beta_2 

DOUBLE PRECISION alpha_iter, QL2_iter 

c Define Common Blocks 


COMMON /MODEL2/ alpha, gamma, lambda2, lambda3. Cl, C2, C3, 

$ Hbf, Hv_2 , Beta_2, QL2, 

$ Ts_2 , Tve_2 

COMMON /GENERAL/ sigma, epsilon, hfg, Cp_liquid, Cp_vapor, Tboil, 
$ Tinitial, Tbf, Te_r, Te_h, hbar_e, 

$ kw, L, hbar__v, P, Aw 

COMMON / ITER_VAL/ EPS_1, xL_l , xR_l , EPS_2 , xL_2 , xR_2 , EPS_3 , 

$ xL_3 , xR_3 
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INPUT . DAT 


$MISC 


sigma 

= 

5.67E-8 

epsilon 

= 

1 . DO 

hfg 

= 

2257. E3 

Cp_liquid 

= 

4217 .DO 

Cp„vapor 

= 

2060 . DO 

Tboil 


373. DO 

Tinitial 

= 

298. DO 

Tbf 

= 

393. DO 

Te_r 

= 

1148. DO 

Te_h 

= 

1148. DO 

hbar_e 

= 

1000 


$ 


$MATERIAL 

MatName = 'Aluminum' 

kw = 155. 8D0 

L = 0.012 

$ 

$MODEL2 
hbar_v 
P 

Aw 

$ 


$ITER__PAR 


EPS„1 

= 

ID -5 

xL_l 

r= 

1 . 11111111D0 

xR_l 

- 

9 . 99999999D0 

EPS_2 

= 

ID - 5 

xh_2 

= 

0 . 00111111DO 

xR_2 

= 

0 . 00999999D0 

EPS_3 

= 

ID -5 

xL__3 

= 

0 . 00111111D0 

xR_3 

= 

0.00999999D0 


$ 


_ONLY 
= 300. DO 
- 0.012 
= 0.000036D0 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Define the Variables 


sigma 

epsilon 

hfg 

Cp_liquid 

Cp„vapor 

Tboil 

Tinitial 

Tbf 

Te__r 

Te_h 

hbar_e 

MatName 

kw 

L 


the Stef an-Boltzman Constant [W/mA2*KM] 

emissivity of the surface 

latent heat of vaporization [J/kg] 

specific heat at constant pressure of the coolant liquid 
[ J/kg*K] 

specific heat at constant pressure of the coolant vapor 
[ J/kg*K] 

boiling temperature of the coolant fluid [K] 

initial temperature of the coolant fluid [K] 

backface temperature of the wall [K] 

radiative temperature of environment [K] 

equivalent temperature of environment due to convection 

[K] 

convection heat transfer coefficient of environment 
[W/m A 2*K] 

text description of material [characters <- 50] 
thermal conductivity of wall material [W/m*K] 
thickness of wall [m] 
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c hbar_ 
c P 
C Aw 

c EPS 1 

c 

c xL_l 
c 

c xR___l 
c 

c EPS_2 
c 

c xL_2 
c 

c xR__2 
c 

c EPS__3 
c 

c xL_3 
c 

c xR__3 
c 


convection heat transfer coefficient of vapor [W/m A 2*K] 

- perimeter of one channel [m] 

- area of one node of channel and surrounding wall [m A 2] 
convergence criterion supporting surface temperature 
iteration in Model I 

initial iteration guess supporting surface temperature 
in Model I 

initial iteration guess supporting surface temperature 
in Model I 

- convergence criterion supporting alpha iteration in 
Model II 

initial iteration guess supporting alpha iteration in 
Model II 

initial iteration guess supporting alpha iteration in 
Model II 

- convergence criterion supporting QL2 iteration in 
Model II 

initial iteration guess supporting QL2 iteration in 
Model II 

- initial iteration guess supporting QL2 iteration in 
Model II 
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Example of OUTPUT. OUT 


**************************** 

************************* 

* Aluminum 

* 

* 

* kw = 155.8 W/m*K 

* 

* 

* L = .0120 m 

* 

* Te_r = 1148.0 K 

* 

* Te_h = 1148.0 K 

* 

* Tbf = 393.0 K 

* 

* hbar_e = 1000.0 W/mA2*K 

* 

* hbar_v = 300.0 W/mA2*K 

* 

* P = .0120 m 

* 

* Aw = .00003600 mA2 

* 

* MODEL I 

* 


* 

Ts I 

= 

453.9 K 


* 

* 

mdot/A I 

= 

.3022 

kg/mA2*s 

* 

* 

MODEL II 



* 

* 

Ts II 

= 

453.0 K 


* 

* 

Tve_II 

= 

426.7 K 


* 

* 

mdot/A II 

= 

.2948 

kg/mA2*s 

* 

* 

MODEL III 



* 

* 

Ts_III 

= 

393.0 K 


* 

* 

mdot/A_IIl 


.3311 

kg/mA2*s 

* 

* * 

*********************** 

**************************** 
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